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Perturbation theory predictions and Monte Carlo simulations for the 2-d 0{n) 

non-linear a-models 
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Dipartimento di Fisica dell’Universitd and INFN, Piazza Torricelli 2, 56126-Pisa, Italy 

By using the results of a high-statistics (C9(10^) measurements) Monte Carlo simulation we test 
several predictions of perturbation theory on the 0(n) non-linear cr-model in 2 dimensions. We study 
the 0(3) and 0(8) models on large enough lattices to have a good control on finite-size effects. The 
magnetic susceptibility and three different definitions of the correlation length are measured. We 
check our results with large-n expansions as well as with standard formulae for asymptotic freedom 
up to 4 loops in the standard and effective schemes. For this purpose the weak coupling expansions 
of the energy up to 4 loops for the standard action and up to 3 loops for the Symanzik action are 
calculated. For the 0(3) model we have used two different effective schemes and checked that they 
lead to compatible results. A great improvement in the results is obtained by using the effective 
scheme based on the energy at 3 and 4 loops. We find that the 0(8) model follows very nicely 
(within few per mille) the perturbative predictions. For the 0(3) model an acceptable agreement 
(within few per cent) is found. 


I. INTRODUCTION 

According to perturbation theory (PT), the 0(ri) non-linear ct- model in 2 dimensions for n > 3 resembles Yang- 
Mills theories in 4 dimensions. Both are asymptotically free |Q,|j and present a spontaneous generation of mass. 
Moreover for n=3 the model has a non-trivial topological content |^. Consequently these models are considered as 
good toy models for testing methods and solutions in 4-dimensional Yang-Mills theories. In condensed matter physics 
these models have applications in the study of ferromagnetic systems. 

There is an extensive literature devoted to investigate the validity of PT in these models on the lattice and in 
particular the onset of scaling (see for instance [^-|[). In [||j|] the 0(3) model was analyzed by using improved 
actions. The obtained results still differ from the exact calculated mass-gap by ^ 15%. In ||] the authors made 
use of the perturbative f3 function up to 3 loops |l^. In faster updating algorithms were used. The mass-gap for 
the 0(3) model was calculated in by using the standard action and an overrelaxed alrarithm. Up to a correlation 
length ^ 300 (in units of lattice spacings) it showed a deviation from the exact result p] of about 20%. The 0(4) 
and 0(8) models with standard action were studied in Q by using the cluster algorithm p^. The deviation from the 
exact result for the 0(8) model at correlation lengths ~ 30 was a few per cent. In an analysis of the performance 
of different lattice geometries for the standard action of the 0(3) model was presented. There was no clear signal of 
an earlier onset of asymptotic scaling. 

The use of PT for such models is not guaranteed. The Mermin-Wagner theorem [|^ states that continuous 
symmetries in 2-dimensional theories cannot be spontaneously broken. Therefore PT, which is an expansion around 
a trivial vacuum, is not a priori well-founded. Motivated by this observation and by the lack of clear asymptotic 
scaling in the previous literature, it has been argued that all 0{n > 2) models undergo a Kosterlitz-Thouless 
(ATT) phase transition at some finite beta Pkt- 

In the present work we have performed a high-statistics simulation (O(IO^) measurements) for the 0(3) and 0(8) 
models on the lattice up to correlations ~ 130 for the 0(^ model and ~ 70 for the 0(8) model. For the 0(3) model 
we have used the tree-level improved Symanzik action ||^ and for the 0(8) model the standard action. We have 
measured the magnetic susceptibility and three different definitions of the correlation length and compared the results 
with both the PT and KT set of predictions. We have computed also some scaling ratios which are particularly 
sensitive to the PT versus KT scenarios. We have made use of the corrections to asymptotic scaling in PT up to 4 
loops in both the standard and effective schemes @ for the 0(8) model and up to 3 loops for the 0(3) model. An 
effective scheme can be defined by using any short distance dominated operator; we have used the density of energy 
operator fl^ . Hereafter we will call it indistinctly effective or energy scheme. To include the analysis in this energy 
scheme, new analytic results are reported in this paper: the 4-loop coefficient in the weak coupling expansion of the 
energy for the standard action and the complete calculation of all coefficients up to 3 loops for the Symanzik action. 

We have used two different definitions of energy operators for the Symanzik action and checked that the corre¬ 
sponding effective schemes agree. Lacking a rigorous treatment for these schemes, this check becomes an important 
test. 

We have avoided strong coupling effects by starting the simulations at large enough correlation lengths. The minimal 
correlation was ~ 10 for the standard action and ~ 16 for the tree-level Symanzik action. 
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We have not made use of finite-size scaling (questioned due to the validity of PT whenever the limit p = L/^ ^0 
holds, where L is the lattice size and ^ any characteristic correlation length) and we have used rather large lattices 
(p ~ 7 — 10) in order to control the finite-size effects. We have checked that the finite-size effects at these p values 
are negligible. 

We are able also to compare the large-n predictions with our data. In particular, we have checked the relationship 
between the two correlation lengths and (see eq. (^) below) known up to 0{l/n) and the prediction for 
the magnetic susceptibility, known up to 0{l/n'^). 

In section 2 we will show the predictions of both PT and KT for the model as well as some necessary I/n expansions. 
In section 3 we will describe our simulations and give the results while in section 4 we will compare them with the 
two different scenarios described in section 2. In this section we will also use the Monte Carlo data of ref. for the 
0(3) model with standard action to check the presently known 4-loop perturbative computations. Our conclusions 
are given in section 5. In the appendix we will show some technical details concerning the perturbative computation 
of the energy up to 4 loops for the standard action and up to 3 loops for the Symanzik action. 


II. PREDICTED SCENARIOS FOR THE cr-MODELS 


The 0(n} non-linear cr-model in 2 dimensions is defined formally in the continuum by the action 

S = ^ J (2-1) 

where 4>{x) is an n-component real scalar field, together with the constraint (pix)"^ = 1 for all x. (3 is the inverse of 
the bare coupling constant. On the lattice one can regularize this theory by making use of different actions. For our 
simulation we chose the standard action 

^standard = _/? ^ ■ ${x + fl) (2.2) 

and the tree-level improved Symanzik action 0 

^Symanzik = f^^x) ■ ${x + p) - ^${x) ■ ${x + 2^)) . (2.3) 

We have measured the magnetic susceptibility x defined as the zero momentum correlation function, 

X=^G{xi,X 2), G{xi,X2) = {${0,Q)-(p{xi,X2)) (2.4) 

Xi,X2 


where we have assumed a symmetric lattice of size L with periodic boundary conditions in both directions and called 
xi and X 2 the two coordinates of the point x. We will need also T defined as the correlation function at the smallest 
lattice non-zero momentum 2t: j 

= E ^^”^^'^G{x^,X2) + \ Y. e""“=/^G(xi,a:2)] . (2.5) 

\ Xi,X2 X\,X2 / 

We have made use also of the wall-wall correlation function defined as 

G{xi) = ^YG{xi,X2)- ( 2 . 6 ) 


We have considered three definitions of correlation lengths, the exponent ial one and the second momenta of 
the correlation function and They are defined as (|a;| = \/x\+ x\) 


^exp 

^(2) 

^'(2) 


|x|^oo lnG(a;i,a;2)’ 
2sin7r/L ’ 


1 Y! \x\^G{xi,X2) 
4 Y.G{xi,X2) 


(2.7) 
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where Yl' indicates that the sum runs over —L/2 + 1 < cci,a ;2 < L/2. The operative definition of on a finite 
lattice was the solution of the equation 


G(ti) cosh ( {t 2 - L/2) /f’'p) = G{t 2 ) cosh ( {h - L/2) /r^^), 


( 2 . 8 ) 


for big enough ti and t 2 where G{t) is the wall-wall correlation and t 2 — ti = At with At = 1, 2. As a function of 
ti, the solution of the previous equation displays a long stable plateau for ^®’'p ^ ti ^ 3 ^®’^p. Anyhow, we chose the 
value and error for self-consistently at ti « 2^®’'p. The result is independent of At (both for the Symanzik act ion 
and the standard one) and we selected the value At = 1. In Figure 1 we show an example of solution of eq. (2.8) as 
a function of ti; the plateau is apparent. 

On the other hand, the value for the definition was extracted from the wall-wall correlation function 


^'( 2 ) ^ 


IlY't'^Git) 
2 YG{t) ■ 


(2.9) 


In the large-L limit and coincide. For finite L the three definitions show rather different finite-size behaviours 





FIG. 1. Solution of eq. (^) as a function of ti for At = 1 for the 0(8) model at /3 = 4.8 and L — 120. The parameter ti 
is given in units of lattice spacings (external labels of the horizontal axis) and units of correlation length (internal labels). 


The scaling of these quantities as predicted in perturbation theory in the large-L limit is 


f n I n-2\n-2 


27r/3 ^ 
n — 2, 




k=l 
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( 2 . 10 ) 


X = C^x 


n — 2 
2 tt(3 


n + 1 


exp 



The form of the expression for ^ is valid for all three definitions. The coefficients of non-universal scaling at and bk 
are action-dependent. They are known up to 3 loops for the Symanzik action |n| and up to 4 loops for the standard 
action The constant is definition- and action-dependent (its dependence on the action is exactly calculable 
in perturbation theory up to an universal constant). is known exactly for the exponential definition. With the 
standard action it is P,p^ 

(1)^ 


C^exp — 




r 1-f 


1 

n — 2 


2 - 5/2 


exp - 


2(n- 2) 


( 2 . 11 ) 


The corresponding constant in the tree-level Symanzik action is easily obtained from ( ^.ll| ) by using the exact 
perturbative result 


-^Symanzik 

.^standard 


= exp 


0.2964 n - 0.0920 
n — 2 


The other constants are not exactly known. For the correlation lengths in the large-n limit we have 


C^f(2) — C^{2) — C^exp (1 — 


0.0032 


O 


In the same limit the value of with the standard action is 


Cx = 0.196 I 1 - 


4.267 


( 2 . 12 ) 


(2.13) 


(2.14) 


In ref. [Q there are numerical results for up to Oil/n^), 


= 0.0127 
Cx = 0.103 


(2.15) 


0(3) (standard action) - 
0(8) (standard action) - 

By us ing eq. ( |2.12| ) the value of for the Symanzik action can be obtained; for n = 3 it is = 0.0625. From eq. 
(2.10) we conclude that the ratio 


UPT = 


27r/3 
n — 2 


n — 1 

n — 2 


1 + 0 


(2.16) 


tends to the constant C^jC'^ as the continuum limit /3 —> cx) is approached. The parentheses contain the corrections 
which are known up to 4 loops for the standard action and 3 loops for the Symanzik action. Hereafter we will call 
this ratio the PT ratio. 

The perturbative expansions of the energy for both the standard action and the Symanzik action are calculated in 
the appendix. 

The correlation length for the 0(2) model, when r = Pkt — /3 is positive and small, scales as 

B 


^ = Hexp 


rl/2 


with A and B positive constants. On the other hand the ratio 

X 


Rkt = 


^2-ri 


(2.17) 


(2.18) 


should be constant as we approach Pkt from below. Here r] = 1/4 is the critical exponent. Following the renormal¬ 
ization group considerations of ref. one can show that r = 1/16 [^. Recent numerical analyses fo r the 0(2) 

model p6|-^ have yielded several values for r which are all consistent with the bound Irj ^ 0.1. Eqs. (2.17-2.18) 
are the expected beh aviou r (and consistent with Monte Carlo simulations) for the 0(2) model. From now on we will 
call the ratio in eq. (2.18) the KT ratio. The KT scenario for the 0(n) model is the extension of this behaviour for 
n > 3. 

In ref. a fit of Monte Carlo data for x and ^ for the 0(3) model with standard action was performed. Within 
errors it gave a constant for the KT ratio and a strong decrease far from constant for the PT ratio. We have 
simulated the 0(3) model with the tree-level Symanzik action in order to check the results of ref. with an 
action classically closer to the continuum limit. 
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III. THE MONTE CARLO SIMULATION 


We have performed an extensive Monte Carlo simulation of the 0(3) model with Symanzik action and the 0(8) 
model with standard action. In tables I and II we show the corresponding sets of raw data. The statistical error of 
the three entries for in table II for L = 50,100, 200 and the two entries for L = 150, 300 in table I display a 

strong dependence on the lattice size. This can be explained by taking into account that the definition of, for example 
involves the quantity ^ \x\'^G{x). For a big enough lattice the “signal” is independent on size, while the “noise” 
grows as the volume. A similar argument can be used for 


Table I 

Raw Monte Carlo data for 0(3) with Symanzik action. The second row for L = 300 was used only for checks of 
finite-size dependenee. 


f3 

L 

10 ® • stat 

X 

^exp 



1.40 

150 

8.1 

361.41(26) 

16.216(23) 

16.195(17) 

15.441(17) 

1.40 

300 

8.08 

360.99(27) 

16.168(27) 

16.161(71) 

16.050(52) 

1.45 

200 

8 

587.18(44) 

21.519(30) 

21.443(23) 

20.462(22) 

1.50 

260 

6.24 

972.78(83) 

28.793(50) 

28.730(35) 

27.273(34) 

1.55 

340 

18 

1634.41(83) 

38.668(61) 

38.636(36) 

36.410(41) 

1.60 

450 

2.88 

2777.3(3.6) 

52.72(38) 

52.42(10) 

49.18(25) 

1.65 

600 

4 

4743.1(5.2) 

71.08(34) 

70.93(13) 

66.27(20) 

1.70 

800 

7.25 

8125.7(6.7) 

96.67(37) 

96.28(13) 

89.64(19) 

1.75 

1050 

1.24 

13852.7(27.7) 

- 

130.25(36) 

- 


Table II 

Raw Monte Carlo data for 0(8) with standard action. The first and third rows (L = 50 and L = 200 respectively) 
were used only for checks of finite-size dependence. 


/? 

L 

10 ^•stat 

X 

^exp 



4.6 

50 

16 

145.501(77) 

9.7787(75) 

9.7541(40) 

7.5678(69) 

4.6 

100 

3.76 

149.00(17) 

9.881(14) 

9.864(13) 

9.533(20) 

4.6 

200 

16 

149.029(83) 

9.8624(74) 

9.860(35) 

9.842(14) 

4.7 

no 

16 

177.86(10) 

10.9156(80) 

10.913(12) 

10.543(13) 

4.8 

120 

16 

212.13(12) 

12.0745(88) 

12.071(13) 

11.635(15) 

4.9 

140 

16 

253.11(14) 

13.370(10) 

13.358(17) 

13.005(17) 

5.0 

160 

40 

302.600(85) 

14.806(10) 

14.788(11) 

14.431(10) 

5.4 

220 

40 

620.78(18) 

22.193(28) 

22.199(15) 

21.381(15) 

5.8 

340 

65 

1289.32(30) 

33.526(47) 

33.411(16) 

32.357(21) 

6.0 

290 

3.2 

1854.6(2.5) 

40.933(81) 

40.879(78) 

36.530(87) 

6.1 

320 

3.2 

2236.6(3.0) 

45.52(10) 

45.425(87) 

40.49(11) 

6.2 

360 

3.2 

2689.7(3.7) 

50.21(11) 

50.26(10) 

44.99(12) 

6.3 

390 

3.2 

3239.7(4.5) 

55.77(12) 

55.60(11) 

49.35(14) 

6.4 

440 

3.2 

3909.8(5.4) 

61.94(15) 

61.75(12) 

55.06(15) 

6.5 

480 

3.2 

4686.9(6.5) 

68.36(17) 

68.16(14) 

60.64(16) 


We have updated the configurations with the Wolff algorithm El- We verified the absence of autocorrelations 
in the data for the standard action. For the Symanzik action we have used a generalization of this algorithm |m| ] 
which does not completely eliminate the critical slowing down. According to the measured integrated autocorrelation 
time pH l, we have performed 4 decorrelating updatings for this action between successive measurements. Once these 
decorrelating updatings were done, we explicitly checked the absence of autocorrelations in the data for the Symanzik 
action. We have measured x the three definitions of ^ shown in the previous section. The necessary two-point 
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correlation functions were evaluated by using an improved estimator . 


Table III 

Integrated autocorrelation times r ™2 for the energies iff 2 and size (C^) of the average Fortuin-Kasteleyn cluster 
for a Symanzik-improved action simulated on a lattice size L = 100. 


p 

Ef 


E^ 

_int 

^2 

(C#) 

2. 

0.98252(5) 

23.6(1.6) 

0.7788(2) 

26.4(1.8) 

3464.9 

5. 

1.14782(2) 

29.4(2.2) 

0.91534(4) 

31.4(2.4) 

6148.7 

10. 

1.19949(2) 

40.4(3.5) 

0.95812(4) 

42.7(3.8) 

7321.1 


We have also done a few runs at small physical volumes, p = L/^<Cl, to calculate the energy for the Symanzik- 
improved action at large /3, (see appendix). We have realised that the performance of the extension of the Wolff 
algorithm for Symanzik actions is less effective in this regime. In Table III we give the integrated autocorrelation 
times for the calculation of the energies iff and iff respectively on a lattice of size L = 100 after 7 10® 
measurements for several /3. The integrated autocorrelation times must be compared with the values r™* ~ 4 found 
when p 1 [Q. In table III we also give the size (C^) of the average Fortuin-Kasteleyn cluster [^,|^. At very small 
physical volumes the result of a single Wolff updating is an almost global flip of the entire lattice, thus becoming an 
approximate reflection symmetry of the whole system. From Table III we see that the average cluster size becomes 
larger as j3 increases (the total number of sites in our lattice is 10000). The worsening of the performance of the 
algorithm in the p ^ 1 regime can be traced back to this fact. Such behaviour is also visible if the standard action 
Wolff algorithm is used. 

We have run our simulations at very high statistics obtaining rather small statistical errors. Therefore the systematic 
errors can become relevant and they require a careful analysis. We consider three sources of such errors: the finite-size 
effects, the different constants in front of the scaling for the correlation length and the non-universal corrections to 
asymptotic scaling. 

All observables, (other than the energy at very high /3), have been measured at values of the ratio p ^ 7. For 
the 0(8) model and (3 < 6.0 this ratio was p > 10. With these p values the finite-size effects are few parts per 
mille and we will not consider them. We have checked this assertion by performing a few runs at different values of 
the previous ratio. For the 0(3) model with Symanzik action at /3 = 1.40 we have used the lattice size L = 150, 
300 (p = 9, 18 respectively) as shown in Table I. The values obtained for %, and are compatible for 
both sizes. Only shows a clear size dependence. We have imposed the predicted L dependence |20| obtaining 
^( 2 )(2^) = 3.9(10.6)/p^ and ^'^^^(L) = ^'(^)(oo) — 6.2(7) p exp(—p/2). We see that although the size 

dependence of has an exponential fall-off the presence of the multiplicative p and the large coefficient in 
front of the exponential makes our data for at p^l not reliable enough. Instead the data for f}^-) are good in 
spite of the presence of the power-law 1/p^. The size dependence of the data for is as gentle as for the data. 

On the other hand for the 0(8) model with standard action we have simulated the value /3 = 4.6 at three lattice 
sizes: 50, 100 and 200 (p = 5, 10 and 20). Again only displays clearly a size dependence. Fitting the data 
to an exponential for and a power-law for obtain = ^*^^^( 00 ) — 3.0(10.0)/p^ and ^'^^^(L) = 

^'(^^( 00 ) — 5.7(5) p exp(—p/2). As before, the L-dependence is sizeable only for due to the large coefficient in 



the correlation is less size-dependent than and on the other hand allows a better error determination than for 
the exponential definition, (to evaluate the error of we also measured the cross correlation between y and F). We 
will correct the non-perturbative constant Qexp, eq. ( |2.11 ), by dividing all data by 0.9979(9) and 0.9989(4) for 0(3) 
and 0(8) respectively. 

The corrections to universal scaling are the largest source of errors and will be discussed in the next section. 
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IV. DISCUSSION OF RESULTS 


In tables I and II we show the raw data for x and the three definitions of In the following analysis we will use the 
values for and we will write instead of . As we said in the previous section we shall neglect the finite-size 
effects and introduce a corrective factor to the prediction (2.11) for Cj. 


u 


SJLP 

u 


1.20 
1.10 
1.00 
0.90 
0.80 
0.70 

1.35 1.45 1.55 1.65 1.75 



P 


FIG. 2. The ratio between non-perturbative constants for the 0(3) model with Symanzik action. Circles (squares, 

triangles) stand for standard scheme (Uf-scheme, Uf-scheme). Open (full) symbols mean 2-loop (3-loop) data. The data for 
the E 2 -scheme have been slightly shifted horizontally to render the figure clearer. 


A. The 0(3) model with Symanzik action 


From the data for and eq. ( 2.1Cl| ) we can compute the constant C^. We shall call such constant obtained 
from the Monte Carlo data. If PT i s cor rect and asymptotic scaling holds, this number should be independent of /3 
and equal to the prediction of eq. (|2.1l| ) for n=3. Therefore the ratio should be 1. In Figure 2 we show 

such ratio as a function of (3 by using eq. ( 2.10 ) at 2-loop and 3-loop for both t he standa rd and energy schemes. 
We used two different energy schemes defined by the operators Ef and , (eqs. ( |8.17 ,B.18) of the appendix). The 
respective (Ie are 


Pei 


w 


SI 

1 


15/12- Fif’ 


Pe2 


W 


S2 

1 


i-Ei' 


(4.1) 


The perturbative expansions of E^ and E 2 are given in the appendix and the Monte Carlo values for both operators 
are listed in Table IV. 
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Table IV 


Measured values of the two operators ( 8.1'/ ) and ( 8.1^) for the 0(3) model with Symanzik action. 


f3 

L 

10 ^ • stat 

iff 


1.40 

300 

4 

0.840997(3) 

0.661762(2) 

1.55 

340 

4 

0.890991(2) 

0.703158(1) 

1.60 

450 

4 

0.904603(1) 

0.714419(1) 

1.65 

600 

4 

0.917106(1) 

0.724757(1) 

1.70 

800 

4 

0.928573(1) 

0.734231(1) 

1.75 

1050 

0.36 

0.93917(2) 

0.74299(1) 

5.0 

100 

7 

1.14782(2) 

0.91534(4) 

10.0 

100 

7 

1.19949(2) 

0.95812(4) 

15.0 

100 

4 

1.21650(3) 

0.97222(2) 


Figure 2 displays an asymptotic approach to unity for increasing /3. The data in the standard scheme (circles) differ 
from unity by ^ 15%. This is in accordance with previous numerical calculations of with the tree-level Symanzik 
action |^. However, the lack of asymptotic scaling in the energy schemes (squares and triangles) amounts only to 
2 — 3% at 3 loops. Notice also that the two energy schemes agree fairly well; this fact supports the reliability of these 
schemes. This agreement improves as (3 increases. In the previous section we saw that the systematic error in the 
corrective factor Cj( 2 ) /C^exp was of the order of 1 per mille which is negligible in Figure 2. 



P 


FIG. 3. The non-perturbative constant for the 0(3) model with Symanzik action. The constant is given in units of 
Astandard. Circles (sqnaiBS, triangles) stand for standard scheme (iff-scheme, iff-scheme). Open (full) symbols mean 2-loop 
(3-loop) data. The data for the iff-scheme at 3-loop has been slightly shifted horizontally to render the hgure clearer. 
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4.65 



4.55 




❖ 



4.45 





4.35 '—^^ 
2.5 3.0 


3.5 



4.5 


5.0 



FIG. 4. The PT ratio (2.16) for the for the 0(3) model with Symanzik action. Open circles (full circles, open diamonds) 
stand for the 1-loop (2-loop, 3-loop) result in the standard scheme. Open squares (open up-triangles, open down-triangles) 
stand for the 1-loop (2-loop, 3-loop) result in the Flf-scheme. Full squares (full up-triangles, full down-triangles) stand for the 
1-loop (2-loop, 3-loop) result in the -scheme. Some data have been slightly shifted horizontally to render the figure clearer. 


In Figure 3 we show the constant computed from our Monte Carlo data at 2 and 3 loops in the standard and 
energy schemes. At present there are no available exact predictions for this constant. The l/ri^ calculation p^ ] 
provides 0.0625 for the tree-level Symanzik action. After rescaling with (^.12) this number becomes 0.0127. From 


the 3-loop data in the energy scheme of Figure 3 one can infer the estimate = 0.0138(2) which differs by ^ 8% 
from the large-n calculation. This result can be compared with the estimate of ref. which is 0.0146(11); we see 
that both agree within errors (notice that at 2 loops our result would be 0.0145(3); this error includes the imprecision 
among the iff and iff data). The estimate of ref. Q has been obtained by using finite-size scaling techniques . 

We observe that the 1/n-expansion up to order 0{l/'n3) converges well, hence we expect a better agreement if 
further corrections were added. Finally, notice that the two energy schemes agree much better at 3-loop than at 
2-loop. 

The results for the PT ratio are reported in Figure 4 up to 3 loops. The data in the standard scheme are far from 
constant although, as is known for the Symanzik-improved actions Q, the slope is less steep than for the standard 
action case 1 2^ . The data for the two energy schemes at 2 and 3-loop agree completely. Moreover these data are flatter 
indicating that scaling has possibly set in. Assuming this onset of scaling, we derive from the data at 2 and 3 loops 
in these effective schemes \n{C^/C^) = 4.54(2). Using the prediction (2.11) = 0.01249, we obtain = 0.0146(2) 

in good agreement with the value inferred from Figure 3. To show the result at 3-loop, the corresponding coefficient 
of the gamma function for the tree-level Symanzik action has been used. This coefficient can be obtained from |^] 
after correcting a misprint in their eq. (25): the (27r)"‘ dividing the last term in Zi must be instead (27r)^. We thank 
M. Falcioni for correspondence about this point, |p3. The 3-loop coefficient is thus 


Symanzik 

72 


(27r)^ 


(-2.01-b 1.65 n-f 0.362 n^) 


(4.2) 
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FIG. 5. The KT ratio (2.18) for the 0(3) model with Symanzik action. 


Now we want to test the KT formulae ( 2.17 ) and ( |2.1§| ). A best fit of the data for ^ to eq. ( 2.17 ) is rather unstable. 
This can be unde rstood as follows: assuming that the KT transition point Pkt does exist and it is far away, we can 
expand eq. (2.17) in powers of (3/Pkt obtaining 




BP 


= A' exp (B'P). 


(4.3) 


This equation shows that actually we are htting the combination B' = B/{2P^^), therefore the best ht cannot yield 
reliable information about the precise value of Pkt- However, the fact tha t the previous analysis within PT gave 
rather acceptable results indicates that the linear approximation in eq. (^|^) is good and indeed Pkt is much larger 
than our working /3’s. 


In Figure 5 the results for the KT ratio (2.18) are shown. By using the previous conclusion about the large value of 
Pkt, we have assumed that inside the narrow interval 1.4 < /3 < 1.75 the factor K in (2.18) is almost constant. As a 
consequence we did not consider it. In ref. this ratio, calculated for the standard action, looked almost constant 
with the critical exponent 77 = 1/4. We emphasize that our data have smaller error-bars and so the interval in the 
vertical axis is almost 7 times finer for our data. This fact allows us to see that our result is clearly not constant. We 
have estimated the probability Q that the data in Figure 5 follow a straight line. Q is obtained from the tail of the 
probability distribution, (we have assumed a gaussian distribution for the point ordinates). We have obtained less 
than Q = 0.01 which means that with probability ~ 99% the data do not follow a straight line. We have repeated the 
same analysis after removing the first two points (one can argue that they are still far from the scaling region of the 
KT transition). In this case Q = 0.09 which still indicates that the data do not lie on a straight line with probability 
91%. If the constancy of this ratio was to be a true physical effect then our data for the Symanzik-improved action 
should stay also constant. 

A similar probability calculation shows that also the 2 and 3-loop data in the energy scheme of Figure 4 do not 
follow a straight line (although the 1-loop data in this scheme is essentially flat). We remark, however, that the 
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effective schemes and the loop corrections have flattened out the data in the PT ratio. In contrast, the increase of 
the resolution in the statistics has revealed that the KT ratio is not as flat as claimed in ref. |^ . 

Our results for the 0(3) model in the standard scheme do not confirm either of the two scenarios. The lack of 
asymptotic scaling agrees with previous works using the same Symanzik improved action, (P. However, in the energy 
schemes these data present asymptotic scaling at 3 loops within 2 — 3% for the correlation length as well as an estimate 
for the magnetic susceptibility that is in reasonable accordance with previous numerical simulations |]3^ and the 1 /n 
expansion. The PT ratio in the energy scheme shows a much flatter behaviour than in the standard scheme. Moreover 
the agreement between the two energy schemes is a reassuring result. 

On the other hand, in the KT scenario, we have seen that the scaling law ( ^.18 ) is badly satisfied. This is in 
constrast with the data of for the standard action. 


B. The 0(3) model with standard action 

In Table V we show the Monte Carlo results for the 0(3) model with standard action taken from ref. and the 
Monte Carlo energy, (see eq. (^ ) of the appendix), from ref. ^]. The correlation length data corresponds to the 
exponential definition in eq. (|2.7D, so there is no correction factor in this case. 


Data for the 0(3) model with standard action. 

Ml- 


Table V 

The X o,i^d data has been taken from B the energy data from 


/? 

L (for X, e"P) 

X 

^exp 

L (for E) 

E 

1.50 

256 

176.4(2) 

11.05(1) 

128 

0.601597(16) 

1.60 

256 

448.4(7) 

19.00(2) 

128 

0.635722(10) 

1.70 

512 

1263.7(3.3) 

34.44(6) 

256 

0.664240(5) 

1.75 

768 

2197.(15.) 

47.2(2) 

256 

0.676629(4) 

1.80 

768 

3823.(21.) 

64.5(5) 

256 

0.687953(3) 

1.85 

768 

6732.(25.) 

88.7(5) 

256 

0.698351(3) 

1.90 

1024 

11867.(62.) 

122.7(1.1) 

256 

0.707952(3) 

1.95 

1024 

20640.(310.) 

164.8(5.3) 

128 

0.716928(9) 


The asymptotic scaling analysis for these data was done up to 3 loops in Q while the test for the KT scenario was 
done in Here we want to make use of our new perturbative results for the energy up to 4 loops, (eq. 


^ J.15D of 

the appendix), and the results of 0 to test asymptotic scaling in the energy scheme for the magnetic susceptibility, 
the correlation length and the PT ratio. The energy scheme is defined as 


Pe = 


Wl 

l-E' 


(4.4) 


In Figure 6 we show the ratio /C^. The lack of asymptotic scaling in the standard schemes is apparent and the 
energy scheme does not improve it as dramatically as for the Symanzik action. We see that the 4-loop correction in 
the energy scheme is almost negligible and as a result the departure from asymptotic scaling at 3-loop observed in ||] 
is still present at 4-loop. The lack of asymptotic scaling in this figure is ^ 10% for the energy scheme and 15 — 20% 
for the standard one. 

Figure 7 displays the non-perturbative constant as computed from the Monte Carlo data. The data in the 
energy scheme converge around = 0.0130(5) while the 1/n^ prediction ||^ is 0.0127 and the result of was 
0.0146(11). The result with our data for the Symanzik action was 0.0138(2)7^ The several Monte Carlo results are 
compatible with each other suggesting that the truncation error of the series at order 1/n^ amounts to ~ 8% when 
n = 3. 
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1.45 1.55 1.65 1.75 1.85 1.95 
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FIG. 6. The ratio between non-perturbative constants for the 0(3) model with standard action. Monte Carlo data 

from Q. Open circles (full circles, open diamonds) correspond to 2-loop (3-loop, 4-loop) in the standard scheme; open squares 
(full squares, open triangles) correspond to the 2-loop (3-loop, 4-loop) in the energy scheme. Some data have been slightly 
shifted horizontally to render the figure clearer. 


0.022 

0.018 

§ ^ 0.015 

U 

0.011 

0.007 

1.45 1.55 1.65 1.75 1.85 1.95 

P 

FIG. 7. The non-perturbative constant for the 0(3) model with standard action. Monte Carlo data from |^. Open 
circles (full circles, open diamonds) correspond to 2-loop (3-loop, 4-loop) in the standard scheme; open squares (full squares, 
open triangles) correspond to the 2-loop (3-loop, 4-loop) in the energy scheme. Some data have been slightly shifted horizontally 
to render the figure clearer. 
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FIG. 8. The PT ratio for the 0(3) model with standard action. Monte Carlo data from Open circles (full circles, open 
diamonds, full diamonds) correspond to 1-loop (further corrections) in the standard scheme; open squares (full squares, open 
triangles, full triangles) correspond to the 1-loop (further corrections) in the energy scheme. Some data have been slightly 
shifted horizontally to render the figure clearer. 

Finally we show the PT ratio up to 4 loops for the standard and energy schemes in Figure 8. The data for the 
standard scheme is clearly not constant as already seen in |^. However again the data in the energy scheme is 
particularly good and stable and allows the determination ln((7^/C|) = 4.57(2) in excellent agreement with the 
previous determination by using our data for the Symanzik action (as it should this ratio is independent of the 
regularization used). 

Our results for the PT ratio and the magnetic susceptibility are 4.57(2) and 0.0130(5) respectively. These results, 
obtained by using the standard action, agree with the previous ones extracted with the Symanzik action. Besides, 
the 0{l/ri^) estimate of is in good accordance with our data. The deviation from and the exact result 

is still of the order 10% even after the inclusion of the 4-loop correction in the energy scheme. 


C. The 0(8) model with standard action 

Our Monte Carlo data for the 0(8) model is shown in Table II. Our data agree with ref. for the two values of (3 
that we have in common. In Table VI we give the energy data taken from p^ |. The perturbative expansion for the 
energy is reported in eq. (|8.16|) of the appendix. 
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FIG. 9. The ratio for the 0(8) model with standard action. Open circles (full circles, open diamonds) stand for 

the 2-loop (3-loop, 4-loop) approximation in the standard scheme. Open squares (full squares, open triangles) stand for the 
2-loop (3-loop, 4-loop) approximation in the energy scheme. 

In Figure 9 we show for the 0(8) model the equivalent of Figure 2. The data converge towards 1 in both the 
standard and energy schemes. The 4-loop energy scheme for the ratio /C^ yields 1 up to 0.5%. 


Table VI 

Energy data for the 0(8) model (from ref. &■ 


p 

L 

E 

4.6 

128 

0.603836(9) 

4.8 

64 

0.620987(10) 

4.9 

128 

0.629018(9) 

5.0 

64 

0.636812(10) 

5.4 

64 

0.664983(9) 

5.8 

256 

0.688885(3) 

6.1 

256 

0.704805(3) 

6.4 

256 

0.719168(2) 


The figure clearly displays that the data approach 1 monotonically as the number of loops increases. An important 
i ssue then is to understand how big the successive corrections are. At leading order in 1/n the coefficients Ofc in eq. 
( |2.10 ) have been computed up to k=8 |pl|] . Comparing with the exactly known coefficients Oi and 02 we see that the 
large-n approximations 01^2 = 8) are correct up to 90% and 60% respectively US- Assuming a corrective factor 
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MC 


/fc = 1 — 2 such that ak = fk ■ o,k(n = 8 ) for all k, then one can see that the next corrections are small and that the 
convergence towards 1 in Figure 9 is meaningful. 

Figure 10 displays the magnetic susceptibility constant as extracted from the Monte Carlo data, We do not 

show the data in the energy scheme at 2 loops as they are very big 0.108) and would expand too much the vertical 
scale of the figure. Data tend to converge around the value « 0.102. Taking the results at 4-loop in the energy 
schemes we obtain = 0.1028(2). The large-n prediction is 0.0915 (up to 0{l/n), [^) and 0.103 (up to 0{l/n^), 
[p^). This Olljn?) estimate agrees with our result within < 0.5% which is the same amount of deviation from unity 
seen in Figure 9 for the correlation length. Therefore the 1/n expansion agrees fairly well with our data. Notice that 
data in the standard scheme do not converge monotonically; indeed we h ave the sequence “2-loop” > “4-loop” > 
“3-loop”. This is due to the fact that the coefficients bi and 62 in eq- (H^) have opposite signs. 


u 


0.1045 

0.1035 

0.1025 

0.1015 

0.1005 

4.5 5.0 5.5 6.0 6.5 

P 



FIG. 10. The non-perturbative constant as extracted from the Monte Carlo data for the 0(8) model with standard 
action. Open circles (full circles, open diamonds) stand for the 2-loop (3-loop, 4-loop) approximation in the standard scheme. 
Full squares (open triangles) stand for the 3-loop (4-loop) approximation in the energy scheme. 
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FIG. 11. The PT ratio for the 0(8) model with standard action. Full circles (open diamonds, full diamonds) stand for the 
2-loop (further loops) approximation in the standard scheme. Open squares (full squares, open triangles, full triangles) stand 
for the 1-loop (further loops) approximation in the energy scheme. The highest order corrections in the energy scheme have 
been slightly shifted horizontally to render the hgure clearer. 


In Figure 11 we show the PT ratio for the 0(8) model. We show this ratio up to 4 loops. We do not show the data 
at 1 loop in the standard scheme because again they lie far from the window shown in the vertical axis. We have also 
omitted the error bars in the further corrections to render the figure clearer. The data stabilize for large enough Inf 
after having included the non-un iversa l corrections. The conve rgenc e is extremely good. The straight horizontal line 
is the prediction (and error) eq. ( 2.16| ) taking the value of eq. (2.11) for and the result 0.1028(2) for from the 
previous figure. Our data gives Rpt = 2.220(5). 
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FIG. 12. The KT ratio for the 0(8) model with standard action. The open circles are the data from our Monte Carlo 
simulation. The solid line is the prediction of perturbation theory at 4 loops. 


Finally Figure 12 shows the KT ratio for the 0(8) model. As in the 0(3) case, we have neglected the variation of 
inside the interval 4.6 < /3 < 6.5. The solid line is the PT prediction for this ratio using eq. (p.lO) up to 4 loops, 


eq. (2.11) for and the result from Figure 10 for C^. We observe that the KT ratio is not constant and that its 


non-constancy is well explained by PT. Notice that the same set of values for and explain well both the PT 
and KT ratios. 

In conclusion, PT works fairly well for the 0(8) model. The data agree with the exact mass-gap ||^ with a precision 
about 0.5%. Analogously the 1/n^ prediction for the magnetic susceptibility is in fair accor danc e with our data within 
the same error. Moreover the PT and KT ratios are well described by the PT formulae (^.1C|). 


V. CONCLUSIONS 


We have done a Monte Carlo simulation for the 0(3) and 0(8) non-linear cr-models in 2 dimensions. The simulation 
was performed with the tree-level Symanzik action for the 0(3) model and the standard action for the 0(8) model. 
We have improved the statistics with respect to previous works and have taken advantage of the recently calculated 
4-loop corrections to scaling [^. We have taken into account the systematic errors coming from the finite lattice 
size and from the different non-perturbative constants for the correlation length In order to reduce them 
we made use of the correlation length data for the definition, eq. (2.7), and calculated numerically the corrective 
factor to pass from Cfex p, eg . (2.11), to C^( 2 ). The result for this numerically calculated factor was in good agreement 
with the 1/n estimate (2.13). 

The ensemble of independent configurations was created with the fast Wolff algorithm, El- The independence of 
the measurements done on these configurations was explicitly verified. However, in the small physical size regime, 
p = L/^ 1 we discovered a worsening in the performance of this algorithm. We argued that this fact can be 

explained by the presence of large Fortuin-Kasteleyn clusters when p <C 1, see Table III. 
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We have also made use of the data of ref. [|j for the 0(3) model with standard action. 

In all cases we tested the perturbation theory predictions in both the standard scheme (expansions in the bare 
coupling 1//3) and the energy scheme (energy modified coupling 1//3_e)- For this purpose in the appendix we have 
computed the weak coupling expansion of the energy up to 4 loops for the standard action and 3 loops for the Symanzik 
action. The fourth loop term in the standard action and the whole expansion up to 3 loops in the Symanzik action 
are new results of the present paper. Moreover, for the Symanzik action, we computed two different operators, Ef 
and E 2 , in order to check the validity of the energy scheme: they should give almost identical results. This check was 
successful (see figures 2-4). 

We saw that the results for the 0(3) model agree fairly well with PT in the energy scheme. The PT ratio leads to an 
almost constant already at 2 loops for both standard and Symanzik actions. This constant was ln(C';^/C'|) = 4.54(2) 
and 4.57(2) for the standard an d Sym anzik actions respectively. The value observed for the non-perturbative constant 
differs from the prediction ( 2.11 ) by almost 2 — 3% in the Symanzik action and > 10% for the standard one. In 
both cases we refer to the results in the energy scheme. Even though these differences are still too large, they are 
much smaller than when obtained from the expansion in the standard 1//3 (~ 20%). The numbers for the constant 
Cx are 0.0138(2) and 0.0130(5) for the Symanzik and standard actions respectively (both in units of Astandard). The 
1/n^ prediction |2j] is 0.0127. Besides, our determinations are in acceptable accordance with the prediction [^]42 


Cx ~ 0.0145. This number has been extracted from the non-perturbative constant Ai obtained in m and the ratio 
between the on-shell and zero-momentum field-renormalization constants q = at large /3. This 

ratio is known in the 1/n expansion to be q = 1 -I- 0.0132/n -!-•••. We have computed this ratio from our Monte 
Carlo data, ^zero-mom ^on-shell j^j^g constant in front of the wall-wall correlation function for 

large separation t 


Git) 


,_,heu exp(£/(2e^P)) 


cosh((t-L/2)/f"P). 


(5.1) 


The value for ^on-sheii pj-gggnted a plateau as a function of t in the interval ^®’'p/2 ^ t ^ 3^®’'p/ 2 and we chose the 
value at t In Table VII we give our numerical result for the ratio q from our data for the 0(3) model with 

Symanzik action as a function of /3. The average is q = 1.0035(18) which is in excellent agreement with the 1/n 
expansion. The fact that q is close to 1 up to few per cent, implies that the estimate « 0.0145 is valid within few 
per cent. We see again a good performance of the 1/n exp ansion even at n = 3. In particular there is a c onsiderable 
improvement from the Oil/n) app roxim ation in eq. (2.14) to the ©(l/n^) order in eq. ( ^.15 ). This fact 
makes us to suspect that also in eq. ( 2.13 ) the ©(l/n^) term would notably improve the agreement with our numerical 
result for that ratio. 


Table VII 

Results for the ratio q as a function of 13. 



1.40 

1.45 

1.50 

1.55 

1.60 

1.65 

1.70 

q 

1.0052(40) 

1.0056(41) 

1.0028(46) 

1.0031(34) 

1.0007(89) 

1.0043(79) 

0.9981(60) 


Recall that the Symanzik action has been designed to reduce lattice artifacts . However this improvement can 
be overwhelmed by the large corrections to asymptotic scaling. The effective schemes can cure this last problem. 
Hence the combination of an improved action together with the use of an effective scheme should provide the best 
results. This may be the reason for the good agreement between the PT predictions and our data from the 0(3) 
model with Symanzik action within the energy schemes. 

Our analysis of the Monte Carlo resul ts for the 0(8) model reveals a satisfactory agreement between the PT 
predictions and the data. The value ( ^.11 ) for 0{ is recovered within 0.5% and the 0(l/n^) prediction for Cx agrees 
within less than 0.5% with our result 0.1028(2), (again there is a remarkable improvement between the Oil/n) and 
0(l/n^) calculations). Analogously the PT ratio tends to stabilize at ln(0^/0|) = 2.220(5); the same prediction 
calculated from the previous value for Ov and the exact Of (2.11) is shown in Figure 11 as an horizontal line at 
InRpT = 2.224(2). 

We have also checked the set of predictions of the KT scenario for the 0(3) model with Symanzik action and the 
0(8) model with standard action. Figures 5 and 12 show the results for these two cases for the KT ratio. None of 
them yield a constant as it happened for the data of ref. . We stress the fact that our data have better resolution 
as the error bars are almost one order of magnitude shorter than in As for the 0(3) model, we showed that 

the probability of having a straight line after eliminating the first two data points in Figure 5 is less than 10%. The 
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situation for the 0(8) model is much clearer: the data are definitively far from constant. In this case perturbation 
theory predicts fairly well the trend of the data, mainly for the largest correlations. It is worth noticing that the two 
ratios, Rkt and Rpx, are well explained with the same set of parameters Oj and Cy- obtained from our analysis. We 
could not draw a similar PT prediction for the Rkt ratio for the 0(3) model like the solid line in Figure 12 because 
the results for and Oj for the 0(3) model had less precision and the KT r atio is rather sensitive to the precision. 

We also tried a fit of the data for the correlation length to the KT law ( ^.17 ). The fit is unstable because the actual 
value for Pkt (if it is finite) is much larger than our working /3’s 

In summary, PT works well if one includes also the non-universal corrections. Only the correlation length data for 
the 0(3) model with standard action still stays far from the (2.11) prediction, although these non-universal corrections 
improve the accordance by a factor of 2. In this respect, we have seen that the energy scheme ||l^ performs very well 
and it is a reliable scheme as explicitly proved by using two different operators with the Symanzik data for 0(3). In 
ref. [Q the authors calculate the non-universal corrections to scaling for the spherical model, discovering that they 
are absent for the energy scheme. We have seen that this good behaviour is almost preserved at low values of n. 
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VII. ELECTRONIC MEMORANDUM 


A progress report of the present work was sent to the Lattice-96 proceedings, |^. While this progress report 
circulated as an hep-lat preprint, a comment appeared which motivated our reply 


VIII. APPENDIX 

In this appendix we sketch the calculation of the energy up to 4 loops for the standard action and 3 loops for the 
tree-level improved Symanzik action. 


A. Standard action 


We define the energy for the standard action as (not summed on fi\) 

E = (0(0) • 0(0 + A)) 

which in the weak coupling expansion can be written as 

1 "^1 ^2 W 3 W 4 

The first two coefhcients wi and W 2 can be straightforwardly computed giving 


Wi = 


(n- 1) 


W2 = 


jn-i) 

32 


( 8 . 1 ) 


( 8 . 2 ) 


(8.3) 
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FIG. 13. Feynman diagrams contributing to the 3-loop coefficient of the free energy. Crosses stand for insertions of the 
measure lagrangian that comes from the Dirac delta in eq. (^.^). 

The order 0{l/p^) coefficient has been computed in (for the 0(3) model it was also calculated in Q and for 
general n in [®). We have checked their result by computing the diagrams for the free energy in Figure 13 and by 
making use of the relationship 




Z = j V$[x) S{${xp - 1) exp _ 


(8.4) 


V is the space-time volume and T> the standard functional measure. In the evaluation of the Feynman diagrams the 
following identity is useful 

{P1+P2) +{Pl+P 3 ) +{Pl+Pi) = P\ + P2 + pI + P\ — 

Ziijkl = ^^PifiPjfiPkfiPlfit (^- 3 ) 


provided that Pi + P 2 + P 3 + Pi = 0 [Q. We make use of the standard notation, p^ = 2 sin(p^/2) and = J^fj, P^i,- 
Another relation useful during the evaluation of tadpole diagrams is 


{P 1 +P 2 ) =pI+pI 


^ Pi P 2 + odd terms, 


valid for any pair of momenta pi and p 2 ■ 
The result for W'i is 


Wz 


- 1 )' jv I (^ - 1 ) 
16 16 


--K+-J 
6 3 


K and J are finite integrals 


( 8 . 6 ) 


(8.7) 
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K = 


J = 



Do 

^ pI pIpIpI 

Do 

-*-^0 ^9 •'9 ^9 -^9 

Pi P 2 Pi Pi 


0.0958876 

0.136620 


( 8 . 8 ) 


where the measure D 3 is 


D 3 


d^pi P?p2 d^p3 d^p4 
(27r)^ (27r)^ (27r)^ (27r)^ 


( 27 r)^ 6 {pi +P2+P3+ Pi) 


and 

Ay = (pi+pj) -pj-p), 
Ai_j = {pi-pj) -f,-p], 

Ai-j will be used later. 


(8.9) 


( 8 . 10 ) 


FIG. 14. 
crosses. 



Feynman diagrams contributing to the 4-loop coefficient of the free energy. Same meaning as in Figure 9 for the 


In Figure 14 we show the diagrams needed for the evaluation of W 4 . Again eqs. 
identities among momenta are needed. The result is 


and (| 8 . 6 | ) are useful. No new 


1 1 1 1 

' 128 2 4 3 24 8 4 
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( 8 . 11 ) 


3{n-lf 

8 





(»- 1)^ 

32 


Hs- 


K and J are given in eq. ( ^ ) while Hi, ..., are genuine 4-loop integrals 


Hi 


H2 

Hz 

Hi 

Hz 



Ai 2 A 34 Si256 
Pi P 2 Pz Pi Pz pi 

A34 S1234 Si256 
^^2 ^^2 ^2 ^^2 ^ ' 

Pi pi pi Pi pi pi 

Ai 3 A 45 A 2_6 

Pi Pi pi Pi pi pi 

S1234 S3456 Si 256 
^2 ^^2 ^^2 ^^2 ^ ^2 

PI P 2 P 3 Pi Pi Pi 

Ai 2 A 34 A 56 
^2 '^2 ^2 ^2 ^^2 '^2 
PI pi pi pi pi pi 


0.0378134 

-0.0322778 

-0.0136824 

= 0.0411085 

-0.0501528 


The measure for the 4-loop integrals is 

n = d^pi d^p2 d^p3 d^p4 d^ps d^pe 
“ (27r)^ (27r)^ (27r)^ (27r)^ (27r)^ (27r)^ ^ 

(27r)^ 5{jpi +P2+PZ+ Pa) (27r)^ 6{pz + Pe + P3 + P4)- 

Numerically at 4 loops the expansion (|8.2|) reads 


E{p) = 1- 


n — 1 n — 1 


4/3 32/32 

0.00291780 (n - 


0.00726994 (n - 1) -f 0.00599298 (n - l)^ 

1) -f 0.00332878 (n - 1)^ + 0.00156728 (n - 1)^ 


For n = 3 and n = 8 the expansion 


(B.14|) becomes 


E{p, n = 3) = 1 - ^ 
E{I3, n = 8) = 1 - ^ 


1 


I 6/32 

7 

3^ 


0.03851 0.03169 

0.3445 0.7211 

“03 


( 8 . 12 ) 


(8.13) 


(8.14) 


(8.15) 

(8.16) 


B. Symanzik action 

As for the Symanzik action, we have used two different local operators to define the so-called energy-scheme (not 
summed over p\) 


Ef = {p{Q) ■ ${0 + fl)- 1^(0) • ${Q + 2/1)), (8.17) 

Ei = {${ 0)-${0 + p)). (8.18) 

The first operator is the energy density for the Symanzik-improv ed a ction, hence its weak coupling expansion can be 
computed by evaluating the free energy and making use of eq. ( ^.4| ). In ref. Q it was computed up to 2 loops for 
the n = 3 case. We have checked their result which for any n can be written as 
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(8.19) 


-I r SI ^ 

^ 12 p ^2 p3 

» f ‘ = IL ^, 


„S1 


(n-i) 

48 




Yi is a 1-loop integral. The notation lip will mean the inverse propagator for the Symanzik action 


lip ~ P 12 




= J 2 pt- 


The 1-loop integral is 


n = 


d^n Dp 

/ —4 = 2.043576 

/-. (27r)" Hp 


The 3-loop coefficient can be obtained by evaluating a set of diagrams analogous to the one in Figure 13. 
identities are 

Ilp+g -l- Ilp+fe -|- IIp+p = lip + Ilg + life -|- lip — S'®, 

^ =2 ^ ' 2pp2(j'^2fcp2rp, 


valid whenever p-\-q + k + r = Q and 
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for any pair of momenta p and q. The result for rcg is 


<,1 (n- 1) o (n- 1) / 1 1 o 1 o 1 
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Y 2 is a 1-loop integral 


The 3-loop integrals are 


f + TT 


Y 2 = 


d^P °p 


{2Trf (Hp)" 


= 
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f+Tr 


D, 


Af.Af4 


iipi np2 iipj iip^ 


f + TT 


D, 




Hpi iip2 iip3 iip4 


= 4.783071 


= 0.0673316 

= 0.104551 


where Afa = (npi+p 2 - IIpi - np 2 ). 
Numerically Ef is 


15 n - 1 0.0244486 (n - 1) 


For n = 3 it is 


El W = 4^ p2 

0.00449054 (n - 1) + 0.00420822 (n - 1)^ 

■ 

15 1 0.04890 0.02581 
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FIG. 15. Feynman diagrams contributing to the 3-loop coefficient of the operator E 2 , eq. (B.18). Crosses and black spots 
stand for insertions of measure and the operator Sf respectively. 


The second operator used is eq. (8.18|). The computation of the coefficients in the weak expansion 


nS 2 


S2 

2 


„S2 




(8.29) 


requires the evaluation of diagrams with an insertion of the operator in (3.18). In Figure 15 we show the diagrams 
necessary for the 3-loop coefficient. The results for all coefficients are 
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The integrals are 


13 = 
^ = 

j5 = 
js = 


d^P al 


i2nf {Upf 


= 11.816615 
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D. 
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^12 
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S1234 

llpi 
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llpg llp^ 



pI 


= 0.0578002 

= 0.0572726 


'-TT (iipi) iipj iipg n. 


= 0.0809553 

= 0.0867806 


(8.31) 
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Numerically is 


El{l3) = 1 - 


0.207425 (n - 1) 0.0189010 {n - 1) 


/3 /32 

0.00353381 (n - 1) + 0.00354656 (n - 1)^ 

■ 


(8.32) 


For n = 3 it is 


E 2 iP, n 


3) = 1- 


0.41485 


0.03780 0.02125 

“52 


(8.33) 


Another method to calculate the previous coefficients has been proposed in ||5l|-|5^ . The Monte Carlo determination 
of any operator at large /3 can be straightforwardly compared to its perturbative expansion, allowing an estimate of 
the perturbative coefficients. In the last three rows of Table IV we give the values of iff and iff for /3 = 5, 10, 15. 

The 0{1/P) coefficient can be obtained comparing the energy at /3 = 15 with the expression 15/12 — wf^/P and 
1 — /p. We obtain wf^ = 0.502(1) and = 0.4167(3). 

Assuming that the exact first order coefficient is known, one can use the value of the energy at /3 = 10 to determine 
the 0{\/P"^) coefficient obtaining = 0.051(2) and = 0.039(4). 

Similarly, by using the exact two first coefficients and the value at /3 = 5 one obtains wf^ = 0.028(2) and = 
0.022(5). 

These results are clearly influenced by the next orders and likely also by the small size {L = 100) of the lattice 
used to calculate the energies for these large /3’s. A better analysis must use a global fit for all coefficients and higher 
precision in the Monte Carlo determination of the operator. Here we have used this technique just as an approximate 
check for our analytical computation. 
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